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Excited-state self-energy effects in the electronic structure 
of Cu, a prototype weakly correlated system containing states 
with different degrees of localization, are investigated with 
emphasis on the unoccupied states up to 40 eV above the 
Fermi level. The analysis employs the experimental quasipar- 
ticle states mapped under full control of the 3-dimensional 
wavevector k using very-low energy electron diffraction for the 
unoccupied states and photoemission for the occupied states. 
The self-energy corrections to the density-functional theory 
show a distinct band- and k-dependence. This is supported 
by quasiparticle GW calculations performed within the frame- 
work of linearized muffin-tin orbitals. Our results suggest 
however that the GW approximation may be less accurate 
in the localized (i-bands of Cu with their short-range charge 
fluctuations. We identify a connection of the self-energy be- 
havior with the spatial localization of the one-electron wave- 
functions in the unit cell and with their behavior in the core 
region. Mechanisms of this connection are discussed based 
on the local-density picture and on the non-local exchange 
interaction with the valence states. 



I. INTRODUCTION 

Many-body exchange-correlation effects in the inhomo- 
geneous interacting electron system of solids, reflected by 
the quasiparticle E(k) band structure observed in the ex- 
periment, are still far from complete understanding. The 
standard Density Functional Theory (DFT) accounts for 
these effects only in the ground state, yielding the static 
properties of the solids such as the charge density or 
the total energy. Of the DFT eigenvalues e(k), how- 
ever, only the highest occupied one yields the correct 
excitation energy by a DFT analogue to the Koopmans 
theorem, 1 while the others, strictly speaking, have no 
physical meaning as energy levels or excitation energies. 
Description of the quasiparticle excited states, created 
by external photon, electron, etc. impact, is by far more 



difficult (see, e.g., [2,3]). Their energies E(k) deviate 
from e(k) owing to the difference of the excited-state dy- 
namic exchange-correlation potential, or the self-energy, 
from the ground-state static exchange-correlation poten- 
tial Vxc used in the DFT. The former is described by 
a complex energy- and k-dependent non-local self-energy 
operator S, whose imaginary part MS describes the life- 
time broadening of the excitation. Whereas for materi- 
als with strong local on-site correlations the excited-state 
many-body effects are usually described by model Hamil- 
tonians, this approach cannot be extended to weakly 
correlated materials with long-range charge fluctuations. 
Only a few perturbational low-order approaches, such as 
the GW approximation, 2 ' 4,5 are suggested for this case. 

Cu has since long been believed to be a prototype 
weakly-correlated metal for which the DFT ground-state 
picture provides an accurate description of the excita- 
tion spectra (see, e.g., the review [6]). Only recently, 
based on comparison of new angle-resolved photoemis- 
sion (PE) data obtained under full control of the 3- 
dimensional wavevector k with state-of-art DFT calcula- 
tions, we have demonstrated that the self-energy correc- 
tions AS = EQs)- e(k) to the DFT in the valence band 
of Cu display a clear k- and band dependence, reaching 
values as large as 0.5 eV in the d-bands. 7 These results 
gave a serious indication for the importance of self-energy 
effects even in such a supposedly simple metal as Cu. 
In fact, the previous misconception had arisen only due 
to errors of the early DFT calculations which acciden- 
tally matched AS. 8 Our preliminary GW calculations 9 
have confirmed the general trends of the observed self- 
energy behavior. The most recent GW calculations 10 
performed within the pseudopotential framework have 
yielded a striking agreement with the experimental va- 
lence band of Cu concerning the k- and band dependence 
of the AS shifts as well as their absolute values. 

While the self-energy effects in the valence band can 
be studied using PE spectroscopy, the unoccupied states 
still remain largely unexplored. Conventional spectro- 
scopies of the upper states using X-ray absorption or 
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Bremsstrahlung phenomena 11 provide only k- integrated 
information such as the matrix element weighted density 
of states. The k-resolving techniques such as inverse PE 
still involve two electron states and suffer from uncertain- 
ties in the surface-perpendicular wavevector k±. Only 
recently it was realized that in the energy region above 
the vacuum level the quasiparticle E(k) can be mapped 
under full control of the 3-dimensional k using Very-Low- 
Energy Electron Diffraction (VLEED) (see [7,12,13] and 
references therein). As the PE final states are the time- 
reversed LEED states, 14 this information can then be 
used in PE spectroscopy to achieve 3-dimensional map- 
ping of the valence band £(k). 7 ' 13 

Here, we extend the self-energy analysis to the unoccu- 
pied states in the energy region from the vacuum level up 
to 40 eV above the Fermi level E F based on the VLEED 
experimental data. Contrary to a commonly accepted 
view, in this energy region the AS shifts also show a 
clear k- and band dependence. We support our findings 
by GW calculations performed within the framework of 
linearized muffin-tin orbitals, which include the valence 
and core states on equal footing. We show that the AS 
anomalies correlate with the spatial character of the one- 
electron wavefunctions </>k, and endeavour to identify the 
essential physics of this. Our results on the valence band 
suggest that the GW approximation underestimates the 
self-energy shifts due to the short-range charge fluctua- 
tions characteristic of the rf-bands of Cu. 



II. EXPERIMENTAL PROCEDURE AND 
RESULTS 

Our analysis employs the results of the recent VLEED 
and PE measurements performed on the Cu(110) 
surface. 7 The information about the bulk band struc- 
ture, reflected in the VLEED spectra, was insensitive to 
a slight oscillatory relaxation of this surface. 15 This was 
confirmed by measurements on different surfaces 7 as well 
as dynamical VLEED calculations. 16 

Mapping of the unoccupied quasiparticle E(k) using 
VLEED (see [7,12,13] and the references therein) is based 
on the fact that the diffraction process is connected with 
E(k±) along the Brillouin zone (BZ) direction deter- 
mined by parallel momentum conservation (external K|| 
= internal k||). In particular, the extrema in the differ- 
ential VLEED spectrum dT/dE of the elastic electron 
transmission reflect the critical points (CPs) in the k± 
dispersions of the bands whose Bloch wave </>k efficiently 
couples to the incident plane wave. Ky dispersion of 
the dT/dE extrema, reflecting the CPs with k± lying on 
a surface-parallel symmetry line of the BZ, yields then 
E(k) along this line resolved in the 3-dimensional k. In 
our case the angle-dependent VLEED data with K|| vary- 
ing in the TY azimuth of the surface BZ yielded E(k) 
along the TX line of the bulk BZ. 7 

The direct mapping of the dT/dE extrema reflects 



however the CPs in the quasiparticle E(k) near the sur- 
face, corresponding to 0k excited by the VLEED beam 
and damped towards the crystal interior due to the fi- 
nite electron lifetime, expressed by the electron absorp- 
tion potential Vi = ImY>. Such E(k) is smoothed 
in k± (the dispersions in k|| remain nevertheless un- 
smoothed by virtue of the surface-parallel invariance of 
the VLEED process) compared to the quasiparticle -E(k) 
in the bulk. 7 ' 12 The CPs appear as the points of extremal 
(inverse) curvature of the the k± dispersions, and are in- 
trinsically shifted from the corresponding CPs in the bulk 
E(k), recovered in the Vi = limit, within a few tenths 
of eV. To improve the accuracy of AS evaluation, based 
on comparison with the DFT implying Vi — 0, we cor- 
rected our VLEED data from the 1^-induced shifts using 
a model calculation. 

The model calculation employed the empirical pseu- 
dopotential method including Vi. E(k) and 0k with com- 
plex kj_ values were obtained by solving the Schrodinger- 
like equation 

(~£ V2 + Vps ~ iV *) ^ = E ^ k 

where V ps is the pseudopotential. The VLEED spectra 
were calculated within the matching approach (see de- 
tails in [12]). In principle, matching calculations directly 
on top of the DFT calculations would be more straight- 
forward, but this would require more complicated tech- 
niques such as the k • p expansion. 17 In addition to the 
standard matching formalism, we included the surface 
barrier as an additional layer with a skewed-cosine like 
potential and performed matching on the vacuum-barrier 
and barrier-crystal planes. Compared to the step-like 
barrier, this improved the VLEED spectral structure en- 
ergies by ~ 0.3 eV. The shifts between the Vi — CPs 
and dT/dE extrema, found in the model calculations, 
were then used to correct the experimental points back to 
the Vi = limit to reflect the bulk E(k). This procedure 
also compensated the slight shifts of the VLEED spec- 
tral structures caused by the surface barrier and overlap 
of spectral structures. 

The obtained unoccupied quasiparticle E(k) is shown 
in Fig.l compared to the DFT and GW calculations (see 
below). The experimental points reveal only the bands 
with significant coupling to vacuum, and others are not 
seen. We omitted the points where strong multiple-band 
hybridization effects made it difficult to identify the CPs 
on the TX line reliably. The lower Vi-point measured 
by VLEED on the Cu(100) surface 18 is also added. 

Here we will also scrutinize our previous results on the 
valence band E(k) from [7]. They are perfectly con- 
sistent with the whole body of previous experimental 
data compiled in [6], but have a superior accuracy by 
virtue of employing a combined VLEED-PE method in 
which VLEED was used to achieve full control of the 3- 
dimensional k. The experimental E(k) is shown in Fig. 2 
also compared with the calculations. The original data 
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from [7] was processed to include only the peaks relevant 
to the bulk states and compensate the lifetime induced 
intrinsic shifts of the PE peaks notable in the bottom 
of the sp-band (~ O.leV). In principle, the accuracy 
of these data can be further improved by correcting for 
possible intrinsic shifts due to the surface and matrix 
element effects, but this would require demanding calcu- 
lations within the one-step PE theory. 14 



III. DFT AND QUASIPARTICLE 
CALCULATIONS 

Our DFT calculations used the Generalized Gradient 
Approximation (GGA) 19 for Vxc- For Cu the results 
were however almost indistinguishable from the LDA. A 
self-consistent full-potential linearized augmented plane 
wave (FLAPW) method implemented in the WIEN97 
code 20 was used, with the basis extended by local orbitals 
to reduce the linearization errors in the extended energy 
region of the unoccupied states. Spin-orbit coupling was 
also included. It should be noted that nowadays all state- 
of-art calculations on Cu agree within <~ 100 mcV (at 
least in the valence band) which implies that the true 
DFT e(k) is achieved on this energy scale. The AS cor- 
rections to the DFT eigenvalues, given by the expecta- 
tion values Re J n 0£ (S — Vxc) <hs.dr' with the integra- 
tion over the unit cell fi, were calculated in the framework 
of the GW approximation. 2,4 ' 5 In this approximation the 
self-energy is given by 2 

S(r,r';u;) = ^- / du'e^' G(r, r': cu + u')W(r, r'; J) 

In practical calculations, the Green function G is approx- 
imated by a non-interacting one: 



G(r,r';u;)=£ 



u - £k 



In our case, {V'kjEk} are the LDA-DFT wavefunctions 
and eigenvalues. The screened interaction W is given by 

W(r,r';w) = J dr"v(r - r'W', r'; w) 

where the dielectric matrix e is calculated within the 
random-phase approximation without employing the 
plasmon-pole approximation. The present GW calcula- 
tions employed one-electron calculations using the LDA 
exchange-correlation and the Linearized Muffin-Tin Or- 
bitals - Atomic Sphere Approximation (LMTO-ASA) 
method. In contrast to the pscudopotential framework, 
this method allows to treat exchange-correlation with the 
core states on the same footing as the valence states. De- 
tails of the computation may be found in [5] . The calcu- 
lated AS values were added to the FLAPW eigenvalues 
whose computational accuracy at higher energies is bet- 
ter. The GW calculations on Cu are technically rather 
demanding due to the presence of the d-states. 



The results of the DFT and GW calculations are also 
shown in Figs.l and 2. 

IV. SELF-ENERGY CORRECTIONS 

The unoccupied quasiparticlc E(k) in Fig.l deviate 
from the DFT results on the whole by ~ 1 eV. How- 
ever, the upper Ai-band displays an anomalously large 
shift by as much as ~ 2.8 eV. Near T the bands are 
shifted also very unevenly. These anomalies, missed in 
the previous studies, 21 have been resolved here for the 
first time by virtue of the single-state and k-resolving 
nature of VLEED. Reported previously, 7 ' 10 ' 22 a peculiar 
renormalization relative to the DFT is also observed in 
the valence band in Fig. 2, where the shifts of the sp-band 
(up to +0.4 eV in its bottom) and the cZ-bands (about 
—0.5 eV) even differ in sign. In the bottom of the d- 
bands the deviations become smaller, which is due to hy- 
bridization of the d- and sp-bands. 23 Through the unoc- 
cupied and occupied states, the experimental shifts from 
the DFT are thus intriguingly band- and k-dependent. 

The deviations can have two sources. Firstly, the static 
Vxc is not known exactly, and any DFT calculation 
has to resort to a reasonable approximation such as the 
GGA. However, the small difference between LDA and 
GGA results and the excellent description of ground-state 
properties suggest that for Cu this problem is insignifi- 
cant. The second, therefore predominating, source is the 
excited-state self-energy effects due to the difference of 
the dynamic exchange-correlation £ from the static one 
Vxc- Our comparison of the experimental quasiparti- 
cle and DFT energies thus directly yields AS. The AS 
values in the X-point, corresponding in our experiment 7 
to K|| = 0, are given in Table I compared to the GW 
calculations. 

Inclusion of the dynamic exchange-correlation within 
the GW approximation are seen to vastly improve de- 
scription of the experimental excitation energies. In the 
unoccupied states the agreement is almost ideal, in par- 
ticular for the anomalous shift in upper Ai-band, with 
the only clear exception in the T21 point. In the valence 
band the agreement is also excellent for the sp-band. In 
the (i-bands the calculation correctly reproduces the sign 
of the AS shifts, although underestimates their magni- 
tude. Interestingly, previous first-principles many-body 
calculations on Cu using formalisms other than GW (see, 
e.g., [24]) failed to reproduce the AS behavior even qual- 
itatively. 

The observed self-energy anomalies in the unoccupied 
bands should be taken into account in LEED surface 
crystallography (see, e.g., [25]), which so far relied on 
a monotonous energy dependence of AS. In our case, 
for example, the observed shifts of the VLEED spec- 
tral structures would be interpreted as due to a surface 
relaxation, but they are only due to the band- and k- 
dependent self-energy behavior. The observed excellent 
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relevance of the GW approximation for the unoccupied 
states suggests to use it in LEED crystallography. 

V. MECHANISMS OF THE SELF-ENERGY 
BEHAVIOR 

Although the observed self-energy behavior is well re- 
produced by the GW calculations, their physical mecha- 
nisms remain obscured by the heavy computational ma- 
chinery. We will now endeavour to identify the possible 
mechanisms, at least on a qualitative level. 

A. Local-density picture 

The observed anomalies in AS can be traced back 
to different spatial localization (SL) of the one-electron 
Bloch waves 0k, i.e. the distribution of their weight 
within the unit cell. Such a mechanism was originally 
suggested in [22]. It assumes the local-density concept, 
i.e. neglects the non-locality of S. 

We start with a discussion of the AS behavior in 
the homogeneous electron gas. The dynamic exchange- 
correlation S as function of momentum k and electron 
density n was calculated for this model system by Hcdin 
and Lundquist. 26 From their data we derived AS as the 
difference of S to its value at the Fermi wavevector Uf, 
which, by the DFT analogue to the Koopmans theorem, 1 
is essentially the static Vxc- The resulting dependencies 
of AS on E/Ef are shown in Fig. 3 (left). The crucial pa- 
rameter to determine their behavior is the electron den- 
sity. 

In real crystals the electron density is inhomogeneous. 
The electron gas plasmon dip in AS is damped (dashed 
line) by averaging over the varying local density n(r) and 
by interband transitions. Most importantly, as shown in 
Fig. 3 (right), the effective electron density - and conse- 
quently AS - experienced by a one-electron wavefunction 
depends now on its spatial localization: if 0k has large 
weight in the core region where n(r) is high, it will expe- 
rience a large-density energy dependence of AS, charac- 
terized by a strong AS repulsion from Ep- If 0k expands 
well into the interstitial region with low n(r), it will ex- 
perience a small-density dependence with its smaller AS 
repulsion from Ep- Such an SL mechanism goes beyond 
the usual wavefunction localization effect (see, e.g., [27]) 
as not only the strength of the localization (measured by 
the bandwidth) matters but also the region of its local- 
ization. 

In the valence band of Cu the SL mechanism has cer- 
tainly some relevance (see also [22]): The rf-bands local- 
ized in the high-density core region experience large ef- 
fective density and, reflecting the AS behavior in Fig. 3, 
shift lower in energy, whereas the sp-band with its charge 
spread out into the low-density interstitial region expe- 
riences small effective density and shifts in the opposite 



direction. A qualitative support of this picture was ob- 
tained by calculating the effective densities (0k| n ( r )|0k) 
(see [9]), which well correlate with the experimental AS 
values both in the d- and sp-bands. Exclusion of the 
core states from the total n(r) did not affect the correla- 
tion. However, the SL mechanism alone can not provide 
any quantitative description of the self-energy behavior 
mainly due to neglect of the non-locality of the £ opera- 
tor. Indeed, the corresponding calculation 22 returned an 
unrealistic plasmon dip of AS in the unoccupied bands 
and, although properly reproducing the sign of the AS 
shifts in the valence band, severely underestimated their 
magnitude. 

In the unoccupied states, the SL mechanism could not 
be reconciled with the experimental self-energy anomaly 
in the upper Ai-band even qualitatively unless we in- 
cluded into n(r) all core states down to the Is state hav- 
ing a binding energy of ~ 9 keV. Any significant effect 
of such a deep level on the dynamic exchange-correlation 
within the local density picture seems unrealistic. This 
hints on an extremely non-local mechanism of interaction 
with the core states, as we discuss below. 

B. GW picture 

We will now analyze a connection between 0k and 
the self-energy effects within the GW framework, taking 
into account the non-locality of the £ operator acting 
as S0k = J n S (r, r') 0k (r') dr' . First, we introduce a 
local quantity AS (r) comprising local contributions to 
AS, obtained by integrating out the non-locality of X as 

AS (r) = 0£ / {E (r, r') - V xc (r) 5 (r - r')} k (r') dr' 

Integration of AS (r) over the unit cell with the weight 
oc r 2 gives the total self-energy correction 

AS = 4-7r/ r 2 AS (r) dr 
Jn 

The calculated |0 k | 2 , AS (r) and r 2 -weighted AS (r) 
for the valence and unoccupied states in the A-point, 
obtained within the GW framework, are shown in Fig. 4. 
We note, firstly, that the total AS is always formed by 
a balance between positive local AS (r) values in the 
core region and negative values in the interstitial region 
(which is seen better in the r 2 * AS (r) curves). Such 
a behavior of AS (r) is in fact general. Indeed, if E 
is assumed to be diagonal in the LDA Bloch basis as 
£ (r, r') = ^ k 0k(r)Sk0 k ( r ')> fr° m the above equation 
we obtain 

AS(r) ~ (S k -Vxc (r))|0 k | 2 

In the core region (small r) the negative Vxc blows up, 
forcing AS (r) to become large and positive in this region. 
Sk, determined by the particular 0k, acts as a constant 
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negative offset. In the interstitial region (large r) Vxc 
decreases and Ek forces AE (r) to become negative. In 
our case the off-diagonal elements of £ are small, and 
such a pattern holds well for all states despite strong 
variations in the character of 0k- 

At first glance, the SL mechanism seems to revive in 
this picture, with an amendment that the figure of merit 
is in fact the weight of 0k in the core region rather than 
that in the high-density region (these two can somewhat 
differ). Indeed, if the weight of 0k shifts into the core 
region, the positive contribution would increase and the 
quasiparticle level would shift higher in energy. How- 
ever, this is much hampered by the Ek offset, which is 
also sensitive to the character of 0k in a highly involved 
manner. 

In view of the limitations of the SL mechanism, another 
mechanism of the observed self-energy behavior can be 
suggested. Fig. 4 shows that the valence and two unoc- 
cupied Xi-states are different from others in that their 
0k has a large 4s component, blowing up at the nucleus 
(r = 0). These - and only these - states experience an 
anomalous positive AE shift on top of the general trend 
that the AE repulsion from Ep graually increases upon 
going away from Ep. 

To extend this observation to other points in the BZ, 
we scrutinized the FLAPW calculations to determine the 
4s projections (partial 4s charges) of 0k inside the atomic 
spheres, which are given by coefficients in the expansion 
= J2i m AimUim (r) Y im (9, 0). The calculated 4s pro- 
jections are shown in Fig. 5. Their comparison with the 
experimental results in Figs.l and 2 demonstrates that 
the anomalous positive AE correspond everywhere to 
the Ai-statcs having large 4s character. Moreover, in 
Fig. 6 we show the experimental AE values as a function 
of energy and the 4s projections for the whole TX line. 
Again, the anomalous positive AE on top of a regular 
energy dependence correlate everywhere with the large 
4s character of 0k, with the only clear exception in the 
r 2 5'-point. 

The observed connection can not be explained within 
the above local picture, because the r 2 -weighting cancels 
the contribution of AE (r) from the nucleus region to 
the total AE, see Fig. 4. Moreover, other states such as 
X51 and X41, which have even more total weight in the 
high-density region than the states with the 4s character, 
can nevertheless not experience anomalous positive AE. 
The observed connection can only manifest a non-local 
exchange-correlation effect, acting through the Ek offset 
in the AE (r) dependence. 

In the following we propose a qualitative explanation 
for the above anomalous behaviour. Let us focus our- 
selves on the states X 5 > , X 3 and X\ which have almost 
the same correlation energy Ep, as can be seen in Table 
2. It implies that differences in the self-energy shifts can 
only be due the differences in the exchange energy Ex 
and Vxc- The trend in 



E x = - ^ / k (r) 0; (r) v(r - r')0 q (r') k (r') drdv' 
q J 

where q spans the occupied states, can be understood by 
analyzing the character of the states. Indeed, the X51- 
state, which has a large Ap character, has the smallest 
Ex among the three states. This is because the occu- 
pied states have very little Ap character, which makes 
the above exchange integral small. The .^-state, on the 
other hand, has a significant 3d character, which gives 
a large exchange with the occupied 3d bands resulting 
in the largest Ex- The Xi-state, the anomalous one, 
has a large 4s character but the occupied 4s valence 
states are plane-wave like so that Ex is somewhere in 
between those of the X 5 /-state and X 3 -state. The trend 
in Vxc seems clear from the charge distribution of the 
states, with the .X^-state being the most localized in- 
side the muffin-tin sphere. The physical picture arising 
from this qualitative analysis is that the AVstate has an 
anomalously large shift because it has a large 4s charac- 
ter, which makes its exchange relatively small, and at the 
same time its charge distribution has a large weight in- 
side the muffin-tin sphere where the exchange-correlation 
potential is deep, which makes Vxc large. These two ef- 
fects result in a large self-energy shift. We note that the 
amount of 4s character alone is probably insufficient to 
make a quantitative connection to the self-energy correc- 
tion. For example, the lower unoccupied Xi-state is less 
anomalous although it has a large 4s character too. This 
can be understood from the fact that it has a significant 
3d character (Fig. 5), which increases the amount of ex- 
change and reduces the self-energy shift. Such an effect 
of the 3d character is consistent with the negative AE 
shifts of the valence 3d bands. Moreover, our qualitative 
picture neglects the exchange-correlation with the core 
levels, while the recent GW calculations 10 suggest that 
the contribution of the 3s and 3p levels is significant. 

The exchange contribution from the valence states, 
connected with the wavefunction character, explains 
therefore the self-energy anomalies on top of the SL 
mechanism. In the unoccupied states of Cu the exchange 
contribution becomes critical. 



VI. IS THE GW APPROXIMATION ACCURATE 
FOR THE LOCALIZED ORBITALS? 

The recent GW calculation by Marini et al, 10 per- 
formed within the pseudopotential framework, have 
demonstrated a good agreement with the experimental 
valence band of Cu both on the delocalized sp-band and 
the localized ci-bands. They conluded that the GW ap- 
proximation, originally proposed to describe long-range 
charge oscillations, well extends to the localized orbitals 
and short-range correlations found in Cu. 

Our GW results are somewhat different however. They 
also demonstrate almost ideal agreement with the ex- 
periment on the delocalized bands such as the valence 
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sp-band, but on the localized d-bands they tend to un- 
derestimate the negative AE shifts (see also Table 1). 
Interestingly, both calculations yield a wrong sign of AS 
in the bottom of the (i-bands where they hybridize with 
the sp-band. Compared to the pscudopotential frame- 
work, our LMTO-ASA calculations have two fundamen- 
tal advantages: 1) use of the true one-electron wavefunc- 
tions rather than pseudowavefunctions; 2) explicit inclu- 
sion of all core states, not only the 3s and 3p states as in 
[10]. Another recent GW calculation based on the full- 
potential LMTO 28 also yielded small AE shifts in the 
d-bands. The deviations from the experiment, remaining 
in the (i-bands, would therefore indicate certain short- 
comings of the GW approximation in the description of 
localized orbitals and short-range correlations. However, 
any unambiguous conclusions on this point require fur- 
ther analysis, because the observed mismatch between 
different GW calculations is comparable with their nu- 
merical accuracy. 

VII. OTHER MATERIALS 

Experimental k-resolved data on the self-energy effects 
in the unoccupied bands of other materials is still scarce. 
In graphite 29 we have found only a small and regular 
AE shift to higher energies (0.1 to 0.5 eV relative to 
the GGA-DFT) which however notably increased above 
a distinct absorption threshold at 35 eV. A similar sit- 
uation has been found in NbSe2- 17 Such an increase is 
not surprising, because the energy dependence of i?eS is 
linked to that of JmS via the Kramers-Kronig relations. 
Interestingly, it was not reproduced by GW calculations. 
In the valence band of graphite 30 we have identified the 
SL mechanism to cause different AE shifts of the a- and 
7r-states, having different overlap with the in-plane ori- 
ented valence electron density. In [31,32] it was specu- 
lated that widening of the valence band in graphite and 
other non-metals could result from electron localization 
in high-density regions, in agreement with the SL mech- 
anism. 

The unoccupied bands of Ni 33 feature AE anomalies 
resembling those of Cu. It seems that they are typi- 
cal for noble metals. The valence band of Ni (see, e.g., 
[34] and references therein) is characterized by a nar- 
rowing of the d-bands (AE shifts towards Ep) which is 
opposite to Cu. This self-energy effect has a different 
origin - strong on-site correlations in the partially filled 
ci-shell. The narrowing of the sp-band is however simi- 
lar to Cu. Interestingly, in [23] the k-dependence of AE 
in the whole valence band was decomposed into constant 
contributions from the sp- and c?-bands whose weight was 
determined only by the band hybridization. Narrowing 
of delocalized bands has also been observed for simple 
metals like Na. 35 



VIII. CONCLUSION 

Excited-state self-energy effects in the unoccupied and 
valence band electronic structure of Cu were investigated. 
The unoccupied quasiparticle bands in the energy region 
up to 40 eV above Ep are for the first time mapped 
under full control of the 3-dimensional k using VLEED. 
They demonstrate, similarly to the valence bands, a sig- 
nificant band- and k-dcpcndcncc of the AE self-energy 
corrections to the DFT. The observed AE behavior in 
the unoccupied and valence bands correlates with the 
spatial localization of the one-electron wavefunctions in 
the unit cell, and with their behavior in the core region 
expressed by the 4s and 3d projections. These effects 
are described, correspondingly, by the spatial localization 
mechanism based essentially on the local-density picture, 
and by the non-local exchange interaction with the va- 
lence states. In the unoccupied bands of Cu the latter is 
particularly important. Our GW quasiparticle calcula- 
tions performed within the LMTO-ASA framework yield 
almost ideal agreement with the experiment on the de- 
localized bands such as the unoccupied and valence sp- 
bands, but on the localized (i-bands they underestimate 
the AE shifts. This may indicate certain limitations of 
the GW approximation applied to the localized orbitals 
and short-range correlations. The observed band- and 
k-dependent self-energy effects in the unoccupied bands 
have an important implication in LEED studies of the 
surface crystallography. 
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TABLE I. The DFT eigenvalues in the V-point and 
the corresponding AS corrections: experiment (AEbxp), 
our GW calculations using LMTO-ASA (AEgw-lmto), 
and GW calculations using the pseudopotential scheme 
(A'Egw-pp) from [10]. Energies are in eV relative Ef- 
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AEexp 


AI^GW-LMTO 


AT.GW- 




20.34 


2.44 


2.51 




x 3 
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0.45 


0.43 






12.98/13.16 


0.77/0.59 


0.68 




Xi 


7.13 


0.54 


0.81 




x 4 , 


1.39 




0.20 




x 5 


-1.55/-1.39 


-0.44/-0.60 


-0.26 


-0.64 


x 2 


-1.68 


-0.50 


-0.47 




x 3 


-4.51 


-0.34 


0.05 


0.10 


Xi 


-4.98 


-0.29 


0.13 


0.16 



TABLE II. The exchange and correlation energies for un- 
occupied bands in the V-point. All energies are in eV. 





Edft-flapw 


S x 


s c 


Kc 


Xj 


20.34 


-10.07 


-5.89 


-19.83 


x 3 
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-14.91 


-5.90 


-22.46 


Xy 
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-6.64 


-5.68 


-14.19 


x 1 


7.13 


-18.10 


-3.45 


-23.74 


X4I 


1.39 


-10.48 


-2.18 


-13.92 



7 



FIG. 1. VLEED quasiparticle unoccupied bands (dots) 
compared to the DFT e(k) (solid lines). The deviations re- 
veal the band- and k-dependent excited-state self-energy ef- 
fects. The bands with sufficient coupling to the vacuum are 
marked by gray shading. The quasiparticle GW calculations 
are shown by open circles. 

FIG. 2. PE quasiparticle valence band compared to the 
DFT e(k) and to the GW calculations similarly to Fig.l. 

FIG. 3. Mechanism of the spatial localization effect: (left): 
AE for the free electron gas. Its behavior depends on the elec- 
tron density, (right): In real crystals a wavefunction concen- 
trated in the core region experiences a high effective density, 
whereas larger weight in the interstitial region (exaggerated, 
in reality this case implies a more homogeneous distribution) 
leads to a lower density. This results in different AS behavior. 

FIG. 4. |0 k | 2 and GW calculated local AE (r) and 
r 2 -weighted AE (r) for the unoccupied and valence states in 
the X-point. Note the change in the r scale at 0.1 a.u. The 
curves are normalized to set their variation to 1. AE (r) is 
positive in the core region and negative in the interstitial re- 
gion. The X\ states with the anomalous positive AE values 
have the 4s character. 

FIG. 5. Calculated 4s and 3d character of (fik, expressed by 
the corresponding projections within the atomic spheres (in 
electrons), and represented by the size of the circles. 

FIG. 6. Experimental AE vs the DFT energy and 4s char- 
acter of cj>k from Fig. 5, represented by the size of the vertical 
bars. The dashed line is a guide for the eye showing a regu- 
lar energy dependence of AE for the states with negligible 4s 
character. The anomalous positive AE values correlate with 
the 4s character of the one-electron states. 
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